Entropic barrier of water permeation through single-file channels.

Facilitated water permeation through narrow biological channels is fundamental for all forms of life. Despite its significance in health and disease as well as for biotechnological applications, the energetics of water permeation are still elusive. Gibbs free energy of activation is composed of an enthalpic and an entropic component. Whereas the enthalpic contribution is readily accessible via temperature dependent water permeability measurements, estimation of the entropic contribution requires information on the temperature dependence of the rate of water permeation. Here, we estimate, by means of accurate activation energy measurements of water permeation through Aquaporin-1 and by determining the accurate single channel permeability, the entropic barrier of water permeation through a narrow biological channel. Thereby the calculated value for [Formula: see text] = 2.01 ± 0.82 J/(mol·K) links the activation energy of 3.75 ± 0.16 kcal/mol with its efficient water conduction rate of ~1010 water molecules/second. This is a first step in understanding the energetic contributions in various biological and artificial channels exhibiting vastly different pore geometries.

Structural basis of NINJ1-mediated plasma membrane rupture in cell death.

Eukaryotic cells can undergo different forms of programmed cell death, many of which culminate in plasma membrane rupture as the defining terminal event1-7. Plasma membrane rupture was long thought to be driven by osmotic pressure, but it has recently been shown to be in many cases an active process, mediated by the protein ninjurin-18 (NINJ1). Here we resolve the structure of NINJ1 and the mechanism by which it ruptures membranes. Super-resolution microscopy reveals that NINJ1 clusters into structurally diverse assemblies in the membranes of dying cells, in particular large, filamentous assemblies with branched morphology. A cryo-electron microscopy structure of NINJ1 filaments shows a tightly packed fence-like array of transmembrane α-helices. Filament directionality and stability is defined by two amphipathic α-helices that interlink adjacent filament subunits. The NINJ1 filament features a hydrophilic side and a hydrophobic side, and molecular dynamics simulations show that it can stably cap membrane edges. The function of the resulting supramolecular arrangement was validated by site-directed mutagenesis. Our data thus suggest that, during lytic cell death, the extracellular α-helices of NINJ1 insert into the plasma membrane to polymerize NINJ1 monomers into amphipathic filaments that rupture the plasma membrane. The membrane protein NINJ1 is therefore an interactive component of the eukaryotic cell membrane that functions as an in-built breaking point in response to activation of cell death.

Not sorcery after all: Roles of multiple charged residues in membrane insertion of gasdermin-A3.

Gasdermins execute programmatory cell death, known as pyroptosis, by forming medium-sized membrane pores. Recently, the molecular structure of those pores as well as the diversity in their shape and size have been revealed by cryoTEM and atomic force microscopy, respectively. Even though a growth of smaller to larger oligomers and reshaping from slits to rings could be documented, the initiation of the gasdermin pore formation remains a mystery. In one hypothesis, gasdermin monomers insert into membranes before associating into oligomeric pores. In the other hypothesis, gasdermin oligomers preassemble on the membrane surface prior to membrane insertion. Here, by studying the behavior of monomeric membrane-inserted gasdermin-A3 (GSDMA3), we unveil that a monomeric gasdermin prefers the membrane-adsorbed over the membrane-inserted state. Our results thus support the hypothesis of oligomers preassembling on the membrane surface before membrane penetration. At the same time, our simulations of small membrane-inserted arcs of GSDMA3 suggest that the inserting oligomer can be small and does not have to comprise a full ring of approximately 26-30 subunits. Moreover, our simulations have revealed an astonishingly large impact of salt-bridge formation and protein surroundings on the transmembrane passage of charged residues, reducing the energetic cost by up to 53% as compared to their free forms. The here observed free energy barrier of mere 5.6 kcal/mol for the membrane insertion of monomeric GSDMA3 explains the surprising ability of gasdermins to spontaneously self-insert into cellular membranes.

Multiple pore lining residues modulate water permeability of GlpF.

The water permeability of aquaporins (AQPs) varies by more than an order of magnitude even though the pore structure, geometry, as well as the channel lining residues are highly conserved. However, channel gating by pH, divalent ions or phosphorylation was only shown for a minority of AQPs. Structural and in silico indications of water flux modulation by flexible side chains of channel lining residues have not been experimentally confirmed yet. Hence, the aquaporin "open state" is still considered to be a continuously open pore with water molecules permeating in a single-file fashion. Using protein mutations outside the selectivity filter in the aqua(glycerol)facilitator GlpF of Escherichia coli we, to the best of our knowledge, for the first time, modulate the position of the highly conserved Arg in the selectivity filter. This in turn enhances or reduces the unitary water permeability of GlpF as shown in silico by molecular dynamics (MD) simulations and in vitro with purified and reconstituted GlpF. This finding suggests that AQP water permeability can indeed be regulated by lipid bilayer asymmetry and the transmembrane potential. Strikingly, our long-term MD simulations reveal that not only the conserved Arg in the selectivity filter, but the position and dynamics of multiple other pore lining residues modulate water passage through GlpF. This finding is expected to trigger a wealth of future investigations on permeability and regulation of AQPs among others with the aim to tune water permeability for biotechnological applications.

The Hidden Intricacies of Aquaporins: Remarkable Details in a Common Structural Scaffold.

Evolution turned aquaporins (AQPs) into the most efficient facilitators of passive water flow through cell membranes at no expense of solute discrimination. In spite of a plethora of solved AQP structures, many structural details remain hidden. Here, by combining extensive sequence- and structural-based analysis of a unique set of 20 non-redundant high-resolution structures and molecular dynamics simulations of four representatives, key aspects of AQP stability, gating, selectivity, pore geometry, and oligomerization, with a potential impact on channel functionality, are identified. The general view of AQPs possessing a continuous open water pore is challenged and it is depicted that AQPs' selectivity is not exclusively shaped by pore-lining residues but also by the relative arrangement of transmembrane helices. Moreover, this analysis reveals that hydrophobic interactions constitute the main determinant of protein thermal stability. Finally, a numbering scheme of the conserved AQP scaffold is established, facilitating direct comparison of, for example, disease-causing mutations and prediction of potential structural consequences. Additionally, the results pave the way for the design of optimized AQP water channels to be utilized in biotechnological applications.

A cholesterol analog stabilizes the human ß2-adrenergic receptor nonlinearly with temperature.

In cell membranes, G protein-coupled receptors (GPCRs) interact with cholesterol, which modulates their assembly, stability, and conformation. Previous studies have shown how cholesterol modulates the structural properties of GPCRs at ambient temperature. Here, we characterized the mechanical, kinetic, and energetic properties of the human ß2-adrenergic receptor (ß2AR) in the presence and absence of the cholesterol analog cholesteryl hemisuccinate (CHS) at room temperature (25°C), at physiological temperature (37°C), and at high temperature (42°C). We found that CHS stabilized various structural regions of ß2AR differentially, which changed nonlinearly with temperature. Thereby, the strongest effects were observed for structural regions that are important for receptor signaling. Moreover, at 37°C, but not at 25° or 42°C, CHS caused ß2AR to increase and stabilize conformational substates to adopt to basal activity. These findings indicate that the nonlinear, temperature-dependent action of CHS in modulating the structural and functional properties of this GPCR is optimized for 37°C.

Gasdermin-A3 pore formation propagates along variable pathways.

Gasdermins are main effectors of pyroptosis, an inflammatory form of cell death. Released by proteolysis, the N-terminal gasdermin domain assembles large oligomers to punch lytic pores into the cell membrane. While the endpoint of this reaction, the fully formed pore, has been well characterized, the assembly and pore-forming mechanisms remain largely unknown. To resolve these mechanisms, we characterize mouse gasdermin-A3 by high-resolution time-lapse atomic force microscopy. We find that gasdermin-A3 oligomers assemble on the membrane surface where they remain attached and mobile. Once inserted into the membrane gasdermin-A3 grows variable oligomeric stoichiometries and shapes, each able to open transmembrane pores. Molecular dynamics simulations resolve how the membrane-inserted amphiphilic ß-hairpins and the structurally adapting hydrophilic head domains stabilize variable oligomeric conformations and open the pore. The results show that without a vertical collapse gasdermin pore formation propagates along a set of multiple parallel but connected reaction pathways to ensure a robust cellular response.

Monitoring the binding and insertion of a single transmembrane protein by an insertase.

Cells employ highly conserved families of insertases and translocases to insert and fold proteins into membranes. How insertases insert and fold membrane proteins is not fully known. To investigate how the bacterial insertase YidC facilitates this process, we here combine single-molecule force spectroscopy and fluorescence spectroscopy approaches, and molecular dynamics simulations. We observe that within 2 ms, the cytoplasmic α-helical hairpin of YidC binds the polypeptide of the membrane protein Pf3 at high conformational variability and kinetic stability. Within 52 ms, YidC strengthens its binding to the substrate and uses the cytoplasmic α-helical hairpin domain and hydrophilic groove to transfer Pf3 to the membrane-inserted, folded state. In this inserted state, Pf3 exposes low conformational variability such as typical for transmembrane α-helical proteins. The presence of YidC homologues in all domains of life gives our mechanistic insight into insertase-mediated membrane protein binding and insertion general relevance for membrane protein biogenesis.

Native-like membrane models of E. coli polar lipid extract shed light on the importance of lipid composition complexity.

BACKGROUND: Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments. RESULTS: To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly. CONCLUSIONS: The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.

Lipids and Phosphorylation Conjointly Modulate Complex Formation of ß2-Adrenergic Receptor and ß-arrestin2.

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the ß2-adrenergic receptor and ß-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor's intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

Lipid Dynamics in Membranes Slowed Down by Transmembrane Proteins.

Lipids and proteins, as essential components of biological cell membranes, exhibit a significant degree of freedom for different kinds of motions including lateral long-range mobility. Due to their interactions, they not only preserve the cellular membrane but also contribute to many important cellular functions as e.g., signal transport or molecular exchange of the cell with its surrounding. Many of these processes take place on a short time (up to some nanoseconds) and length scale (up to some nanometers) which is perfectly accessible by quasielastic neutron scattering (QENS) experiments and molecular dynamics (MD) simulations. In order to probe the influence of a peptide, a transmembrane sequence of the transferrin receptor (TFRC) protein, on the dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) large unilamellar vesicles (LUVs) on a nanosecond time scale, high-resolution QENS experiments and complementary MD simulations have been utilized. By using different scattering contrasts in the experiment (chain-deuterated lipids and protonated lipids, respectively), a model could be developed which allows to examine the lipid and peptide dynamics separately. The experimental results revealed a restricted lipid lateral mobility in the presence of the TFRC transmembrane peptides. Also the apparent self-diffusion coefficient of the lateral movement of the peptide molecules could be determined quantitatively for the probed short-time regime. The findings could be confirmed very precisely by MD simulations. Furthermore, the article presents an estimation for the radius of influence of the peptides on the lipid long-range dynamics which could be determined by consistently combining results from experiment and simulation.

Closely related, yet unique: Distinct homo- and heterodimerization patterns of G protein coupled chemokine receptors and their fine-tuning by cholesterol.

Chemokine receptors, a subclass of G protein coupled receptors (GPCRs), play essential roles in the human immune system, they are involved in cancer metastasis as well as in HIV-infection. A plethora of studies show that homo- and heterodimers or even higher order oligomers of the chemokine receptors CXCR4, CCR5, and CCR2 modulate receptor function. In addition, membrane cholesterol affects chemokine receptor activity. However, structural information about homo- and heterodimers formed by chemokine receptors and their interplay with cholesterol is limited. Here, we report homo- and heterodimer configurations of the chemokine receptors CXCR4, CCR5, and CCR2 at atomistic detail, as obtained from thousands of molecular dynamics simulations. The observed homodimerization patterns were similar for the closely related CC chemokine receptors, yet they differed significantly between the CC receptors and CXCR4. Despite their high sequence identity, cholesterol modulated the CC homodimer interfaces in a subtype-specific manner. Chemokine receptor heterodimers display distinct dimerization patterns for CXCR4/CCR5 and CXCR4/CCR2. Furthermore, associations between CXCR4 and CCR5 reveal an increased cholesterol-sensitivity as compared to CXCR4/CCR2 heterodimerization patterns. This work provides a first comprehensive structural overview over the complex interaction network between chemokine receptors and indicates how heterodimerization and the interaction with the membrane environment diversifies the function of closely related GPCRs.

Structural Model of the mIgM B-Cell Receptor Transmembrane Domain From Self-Association Molecular Dynamics Simulations.

Antigen binding to B-cell antigen receptors (BCRs) followed by signaling initiates the humoral immune response. The signaling is intimately coupled to nanoclustering of BCRs and their sorting to specific membrane domains, a process that is ruled by interactions between the BCR transmembrane domain and lipids. While the structure of the extracellular domains of BCRs has been resolved, little is known about the configuration of the constituting four immunoglobulin domains spanning the membrane. Here, we modeled the structure of the transmembrane (TM) domain of the IgM B-cell receptor using self-assembly coarse-grained molecular dynamics simulations. The obtained quaternary structure was validated against available experimental data and atomistic simulations. The IgM-BCR-TM domain configuration shows a 1:1 stoichiometry between the homodimeric membrane-bound domain of IgM (mIgM) and a Ig-α/Ig-ß heterodimer. The mIgM homodimer is based on an asymmetric association of two mIgM domains. We show that a specific site of the Ig-α/Ig-ß heterodimer is responsible for the association of IgM-BCRs with lipid rafts. Our results further suggest that this site is blocked in small-sized IgM-BCR clusters. The BCR TM structure provides a molecular basis for the previously suggested dissociation activation model of B-cell receptors. Self-assembly molecular dynamics simulations at the coarse-grained scale here proved as a versatile tool in the study of receptor complexes.

Dynamic processes in biological membrane mimics revealed by quasielastic neutron scattering.

Neutron scattering is a powerful tool to study relaxation processes in biological membrane mimics in space and time. Combining different inelastic and quasielastic neutron scattering techniques, a large dynamic range can be covered: from atomic to mesoscopic lengths and from femto- to some hundreds of nanoseconds in time. This allows studies on e.g. the diffusion of lipids, the membrane undulation motions, the dispersion of sound waves in membranes as well as the mutual interactions of membrane constituents such as lipids, proteins, and additives. In particular, neutron scattering provides a quite direct experimental approach to the inter-atomic and inter-molecular potentials on length and time scales which are perfectly accessible by molecular dynamics (MD) simulations. Neutron scattering experiments may thus substantially support the further refinement of biomolecular force fields for MD simulations by supplying structural and dynamical information with high spatial and temporal resolution. In turn, MD simulations support the interpretation of neutron scattering data. The combination of both, neutron scattering experiments and MD simulations, yields an unprecedented insight into the molecular interactions governing the structure and dynamics of biological membranes. This review provides an overview of the molecular dynamics in biological membrane mimics as revealed by neutron scattering. It focuses on the latest findings such as the fundamental molecular mechanism of lateral lipid diffusion as well as the influence of additives and proteins on the short-time dynamics of lipids. Special emphasis is placed on the comparison of recent neutron scattering and MD simulation data with respect to molecular membrane dynamics on the pico- to nanosecond time scale.

Critical Comparison of Biomembrane Force Fields: Protein-Lipid Interactions at the Membrane Interface.

Molecular dynamics (MD) simulations offer the possibility to study biological processes at high spatial and temporal resolution often not reachable by experiments. Corresponding biomolecular force field parameters have been developed for a wide variety of molecules ranging from inorganic ligands and small organic molecules over proteins and lipids to nucleic acids. Force fields have typically been parametrized and validated on thermodynamic observables and structural characteristics of individual compounds, e.g. of soluble proteins or lipid bilayers. Less strictly, due to the added complexity and missing experimental data to compare to, force fields have hardly been tested on the properties of mixed systems, e.g. on protein-lipid systems. Their selection and combination for mixed systems is further complicated by the partially differing parametrization strategies. Additionally, the presence of other compounds in the system may shift the subtle balance of force field parameters. Here, we assessed the protein-lipid interactions as described in the four atomistic force fields GROMOS54a7, CHARMM36 and the two force field combinations Amber14sb/Slipids and Amber14sb/Lipid14. Four observables were compared, focusing on the membrane-water interface: the conservation of the secondary structure of transmembrane proteins, the positioning of transmembrane peptides relative to the lipid bilayer, the insertion depth of side chains of unfolded peptides absorbed at the membrane interface, and the ability to reproduce experimental insertion energies of Wimley-White peptides at the membrane interface. Significant differences between the force fields were observed that affect e.g. membrane insertion depths and tilting of transmembrane peptides.

The Multifaceted Role of SNARE Proteins in Membrane Fusion.

Membrane fusion is a key process in all living organisms that contributes to a variety of biological processes including viral infection, cell fertilization, as well as intracellular transport, and neurotransmitter release. In particular, the various membrane-enclosed compartments in eukaryotic cells need to exchange their contents and communicate across membranes. Efficient and controllable fusion of biological membranes is known to be driven by cooperative action of SNARE proteins, which constitute the central components of the eukaryotic fusion machinery responsible for fusion of synaptic vesicles with the plasma membrane. During exocytosis, vesicle-associated v-SNARE (synaptobrevin) and target cell-associated t-SNAREs (syntaxin and SNAP-25) assemble into a core trans-SNARE complex. This complex plays a versatile role at various stages of exocytosis ranging from the priming to fusion pore formation and expansion, finally resulting in the release or exchange of the vesicle content. This review summarizes current knowledge on the intricate molecular mechanisms underlying exocytosis triggered and catalyzed by SNARE proteins. Particular attention is given to the function of the peptidic SNARE membrane anchors and the role of SNARE-lipid interactions in fusion. Moreover, the regulatory mechanisms by synaptic auxiliary proteins in SNARE-driven membrane fusion are briefly outlined.

A Coiled-Coil Peptide Shaping Lipid Bilayers upon Fusion.

A system based on two designed peptides, namely the cationic peptide K, (KIAALKE)3, and its complementary anionic counterpart called peptide E, (EIAALEK)3, has been used as a minimal model for membrane fusion, inspired by SNARE proteins. Although the fact that docking of separate vesicle populations via the formation of a dimeric E/K coiled-coil complex can be rationalized, the reasons for the peptides promoting fusion of vesicles cannot be fully explained. Therefore it is of significant interest to determine how the peptides aid in overcoming energetic barriers during lipid rearrangements leading to fusion. In this study, investigations of the peptides' interactions with neutral PC/PE/cholesterol membranes by fluorescence spectroscopy show that tryptophan-labeled K∗ binds to the membrane (KK∗ ∼6.2 103 M-1), whereas E∗ remains fully water-solvated. 15N-NMR spectroscopy, depth-dependent fluorescence quenching, CD-spectroscopy experiments, and MD simulations indicate a helix orientation of K∗ parallel to the membrane surface. Solid-state 31P-NMR of oriented lipid membranes was used to study the impact of peptide incorporation on lipid headgroup alignment. The membrane-immersed K∗ is found to locally alter the bilayer curvature, accompanied by a change of headgroup orientation relative to the membrane normal and of the lipid composition in the vicinity of the bound peptide. The NMR results were supported by molecular dynamics simulations, which showed that K reorganizes the membrane composition in its vicinity, induces positive membrane curvature, and enhances the lipid tail protrusion probability. These effects are known to be fusion relevant. The combined results support the hypothesis for a twofold role of K in the mechanism of membrane fusion: 1) to bring opposing membranes into close proximity via coiled-coil formation and 2) to destabilize both membranes thereby promoting fusion.

Dynamic Cholesterol-Conditioned Dimerization of the G Protein Coupled Chemokine Receptor Type 4.

G protein coupled receptors (GPCRs) allow for the transmission of signals across biological membranes. For a number of GPCRs, this signaling was shown to be coupled to prior dimerization of the receptor. The chemokine receptor type 4 (CXCR4) was reported before to form dimers and their functionality was shown to depend on membrane cholesterol. Here, we address the dimerization pattern of CXCR4 in pure phospholipid bilayers and in cholesterol-rich membranes. Using ensembles of molecular dynamics simulations, we show that CXCR4 dimerizes promiscuously in phospholipid membranes. Addition of cholesterol dramatically affects the dimerization pattern: cholesterol binding largely abolishes the preferred dimer motif observed for pure phospholipid bilayers formed mainly by transmembrane helices 1 and 7 (TM1/TM5-7) at the dimer interface. In turn, the symmetric TM3,4/TM3,4 interface is enabled first by intercalating cholesterol molecules. These data provide a molecular basis for the modulation of GPCR activity by its lipid environment.

Binding Characteristics of Sphingosine-1-Phosphate to ApoM hints to Assisted Release Mechanism via the ApoM Calyx-Opening.

Sphingosine-1-phosphate (S1P) is a lysophospholipid mediator carried by the HDL-associated apoM protein in blood, regulating many physiological processes by activating the G protein-coupled S1P receptor in mammals. Despite the solved crystal structure of the apoM-S1P complex, the mechanism of S1P release from apoM as a part of the S1P pathway is unknown. Here, the dynamics of the wild type apoM-S1P complex as well as of mutants were investigated by means of atomistic molecular dynamics simulations. The potential of mean force for S1P unbinding from apoM reflected a large binding strength of more than 60 kJ/mol. This high unbinding free energy for S1P underlines the observed specificity of the physiological effects of S1P as it suggests that the spontaneous release of S1P from apoM is unlikely. Instead, S1P release and thus the control of this bioactive lipid probably requires the tight interaction with other molecules, e.g. with the S1P receptor. Mutations of specific S1P anchoring residues of apoM decreased the energetic barrier by up to 20 kJ/mol. Moreover, the ligand-free apoM protein is shown to adopt a more open upper hydrophilic binding pocket and to result in complete closure of the lower hydrophobic cavity, suggesting a mechanism for adjusting the gate for ligand access.

Exploring the Formation and the Structure of Synaptobrevin Oligomers in a Model Membrane.

SNARE complexes have been shown to act cooperatively to enable the synaptic vesicle fusion in neuronal transmission at millisecond timescale. It has previously been suggested that the oligomerization of SNARE complexes required for cooperative action in fusion is mediated by interactions between transmembrane domains (TMDs). We study the oligomerization of synaptobrevin TMD using ensembles of molecular dynamics (MD) simulations at coarse-grained resolution for both the wild-type (WT) and selected mutants. Trimerization and tetramerization of the sybII WT and mutants displayed distinct kinetics depending both on the rate of dimerization and the availability of alternative binding interfaces. Interestingly, the tetramerization kinetics and propensity for the sybII W89A-W90A mutant was significantly increased as compared with the WT; the tryptophans in WT sybII impose sterical restraints on oligomer packing, thereby maintaining an appropriate plasticity and accessibility of sybII to the binding of its cognate SNARE partners during membrane fusion. Higher-order oligomeric models (ranging from pentamer to octamer), built by incremental addition of peptides to smaller oligomers, revealed substantial stability and high compactness. These larger sybII oligomers may induce membrane deformation, thereby possibly facilitating fast fusion exocytosis.